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Abstract 

Motivated by our recent work [26] , we generalize this work to the interacting non- 
flat case. Therefore in this paper we deal with canonical, phantom and quintom 
models, with the various fields being non-minimally coupled to gravity, within the 
framework of interacting holographic dark energy. We employ the holographic model 
of interacting dark energy to obtain the equation of state for the holographic energy 
density in non-flat (closed) universe enclosed by the event horizon measured from 
the sphere of horizon named L. 
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1 Introduction 



Recent astrophysical data [1] strongly indicate that the universe is accelerating at 
present. Therefore, it is of paramount importance to explain why this is happening. 
Many theories have been proposed recently that try to address this issue, which has 
become the most important problem in cosmology. Although theories that try to modify 
Einstein equations [2] constitute a big part of these attempts, the mainstream explanation 
for this problem, however, is known as theories of dark energy. The most accepted idea 
indicates that a mysterious dominant component, dubbed dark energy, which has negative 
pressure, leads to this cosmic acceleration, though its nature and cosmological origin still 
remain unknown. 

The combined analysis of the cosmological observations suggests that the universe con- 
sists of about 70% dark energy, 30% dust matter (cold dark matter plus baryons), and 
a negligible percentage of radiation. Although little about dark energy is known, we can 
still propose some candidates to describe it. The most natural candidate to account for 
this acceleration is the cosmological constant, by given the problems associated with it 
we turn our attention to dynamical dark energy models. The dynamical nature of dark 
energy, at least at an effective level, can be originated from various fields, such as a canon- 
ical scalar field (quintessence) [3], a phantom field, that is, a scalar field with a negative 
sign of the kinetic term [4], or the combination of quintessence and phantom in a unified 
model named quintom [5]. The advantage of this combined model is that although in 
quintessence the dark energy equation of state parameter remains always greater than — 1 
and in phantom cosmology always smaller than —1, in the quintom scenario it can cross 
-1. 

The general consensus among theorists is that we can not entirely understand the nature 
of dark energy until a complete theory of quantum gravity is established [6]. The dark 
energy problem may be then in essence a problem that belongs to quantum gravity [6]. 
Thus, a complete theory of quantum gravity should explain the properties of dark energy, 
such that the energy density and the equation of state would be determined definitely 
and uniquely. However, in spite of the fact that we still lack a theory of quantum gravity, 
we can make some attempts to probe the nature of dark energy by making use of the 
holographic principle, which is thought to be present in a final theory of quantum gravity. 
In particular, an interesting attempt within the framework of quantum gravity is the holo- 
graphic dark energy (HDE) proposal [7, 8, 9, 10], which is constructed in the light of the 
holographic principle and therefore possesses some significant features of an underlying 
theory of dark energy. The HDE model has been tested and constrained by various as- 
tronomical observations [11, 12] as well as by the Anthropic Principle [13]. Furthermore, 
the HDE model has been extended to include the spatial curvature contribution, i.e. the 
HDE model in non-flat space [14]. 

On the other hand, it is usually assumed that both dark energy and dark energy only 
couple gravitationally. However, given their unknown nature and that the underlying 
symmetry that would set the interaction to zero is still to be discovered, an entirely 
independent behavior between the dark sectors would be very special indeed. Moreover, 
since dark energy gravitates, it must be accreted by massive compact objects such as 
black holes and, in a cosmological context, the energy transfer from dark energy to dark 
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matter may be small but non-vanishing. 

Furthermore, the coupling is not only likely but may even be inevitable [15]. In addi- 
tion, it might explain or at least alleviate the coincidence problem [16, 17]. It was found 
that an appropriate interaction between dark energy and dark matter can influence the 
perturbation dynamics and affect the lowest multipoles of the CMB angular power spec- 
trum [18, 19]. Thus, it could be inferred from the expansion history of the Universe, as 
manifested in several experimental data. Moreover, it was suggested that the dynamical 
equilibrium of collapsed structures such as clusters would be modified due to the coupling 
between dark energy and dark matter [20, 21]. There is no clear consensus on the form 
of the coupling. Most studies on the interaction between dark sectors rely either on the 
assumption of interacting fields from the outset [17, 22], or from phenomenological re- 
quirements [16]. The coupling not only affects the expansion history of the Universe but 
modifies the structure formation scenario as well because we no longer have p m oc a~ 3 . 
The aforesaid interaction has also been considered from a thermodynamical perspective 
[23, 24] and has been shown that the second law of thermodynamics imposes an energy 
transfer from dark energy to dark matter. 

A pressureless dark matter in interaction with holographic dark energy is more than just 
another model to describe an accelerated expansion of the universe. It provides a unifying 
view of different models which are viewed as different realizations of the interacting HDE 
model at the perturbative level [25]. 

In the present paper we extend our recent work [26] to the interacting HDE scenario in 
a non-flat universe. We utilized the horizon's radius L measured from the sphere of the 
horizon as the system's IR cut-off. The organization of our work is as follows: In section 
2 we construct the cosmological scenarios of non-minimally coupled canonical, phantom 
and quintom fields, in the framework of HDE. In section 3 we examine their behavior and 
we discuss their cosmological implications. Finally, in section 4 we summarize our results. 

2 Interacting non-minimally coupled fields of holo- 
graphic dark energy in non-flat universe 

2.1 Canonical field 

We first consider a canonical scalar field with a non-minimal coupling. This case was 
partially investigated in [27], and recently extended in [26]. The action of the universe is 

(1) 

where k 2 is a gravitational constant. In this action we have added a canonical scalar field 
<f), which in non-minimally coupled to the curvature with coupling parameter Lastly, 
the term Cm represents the matter content of the universe and the term x multiplying it 
accounts for the interaction. 

The presence of the non-minimal coupling leads to the effective Newton's constant 

87r<7 e// = (i - c^vt 1 • ( 2 ) 

We proceed now to calculate the equation of state for the HDE density when there 
is an interaction between the HDE and Cold Dark Matter(CDM) with w m = 0. The 
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continuity equations for the HDE and CDM are 



Pa + 3if (1 + w A )p A = -Q, (3) 
p m + 3H Pm = Q (4) 

here H = a/a is the Hubble parameter. The interaction is given by the quantity Q = Tp A 
and describes a decay of the HDE component into CDM with the decay rate given by T. 
Taking the ratio of the two energy densities as r m = p m / p\, the above equations lead to 



fm — 3Hr ri 

Following Ref. [28], we define 



1+7V 

w A H 



3H 



(5) 



Thus, the continuity equations can be written in their standard form as 

p A + 3H(1 + wf)p A = , (7) 

p m + 3H(l + w c ^)p m = . (8) 
We now consider the non-flat Friedmann-Robertson- Walker universe with line element 

rfs 2 = _ dt 2 + fl2(t) ^-^— 2 + r 2 dn 2 ^j . (9) 

where k denotes the curvature of space k — 0,1, —1 for flat, closed and open universe, 
respectively. It must be noticed that a closed universe with a small positive curvature 
(ilk ~ 0.01) is compatible with observations [29, 30]. As usual, we use the Friedmann 
equation to relate the curvature of the universe to the energy density. In the interacting 
case we are dealing with, the Friedmann equations and the evolution equation for the 
scalar field are 

K 2 (p m + PA + \4> 2 + 6^H<M)) 

R2 = o n 7 (10) 

3 (1 - £<pK 2 (j) 2 ) 



+ 3iJ0 + 6^ (H + 2H 2 



-Q (11) 



Pm. + Pa + 3H (p m + pa + Pm + Pa) = 0. (12) 

In these expressions, p m and p m are the pressure and density of the matter content of 
the universe, respectively. Finally, pa and pa are the corresponding components of dark 
energy, which are attributed to the scalar field. 
We define as usual 

Q _Prn_ Pm Q _ PA _ PA Q _ _±_ . , 

m ~ p cr ~ 3M p 2 ^ 2 ' ~ p cr - 3M p 2 F 2 ' h ~ am 2 [ 6) 

where M 2 = (torG eff ) _1 . 

Now we can rewrite the first Friedmann equation as 

il m + il A = 1 + il k . (14) 



1 



This allows us to express r m and = pk/pA in terms of fl A and £\ as 



1 — £l A + i\ ^ _ /^gx 



In the non-flat universe, our choice for HDE density is 
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p A = -(l-^V)£- 2 , (16) 



where L is defined as 

L = ar{t), (17) 
here, a, is the scale factor and r(t) can be obtained from the following equation 

r{t) dr f°°dt R h 



Vl - kr 2 Jt a a 

where Rh is the event horizon. Therefore while Rh is the radial size of the event horizon 
measured in the r direction, L is the radius of the event horizon measured on the sphere 
of the horizon. For the closed universe we have (the same calculation is valid for the open 
universe by a transformation) 

r(t) = —=siny. (19) 
V k 

where y = \/kRh/a. As in [26] we are interested in extracting power-law solutions of our 
cosmological model (10)-(12), in the case of a dark-energy dominated universe (p m ,p m <C 
1). Thus, we are looking for solutions of the form 

a(t) = a t r 

<t>(t) = for*. (20) 

The insertion of the second ansatz allows us to express the HDE density as 

p A (t) = - 2 {l-^K 2 fo 2 t 2s *)L-\ (21) 

Using the definitions Q A = j^- and p cr = 3M 2 H 2 , we get 

HL = (22) 
Now using Eqs.(17, 18, 19, 22), we obtain 

L = HL + ar{t) = - cosy. (23) 

V ^A 

By considering the definition of Eq.(21) for the HDE p A (t), and using Eqs.(22, 23) one 
can find 

p A = -2s 4> r 1 6A ( ^ffZs ) -2H(l- v^cosy) p A . (24) 
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If we substitute the above relation into Eq.(3) and use the definition Q = Tp A , we arrive 
at 



w A (t) 



cos y H 

3 y 3H 



Here as in Ref.[31], we choose the following relation for the decay rate 

T = 3b 2 (l + r m )H 

with the coupling constant b 2 . Using Eq.(15), the above decay rate yields 

r = aWrfi+fW 

"A 



(25) 



(26) 



(27) 



Substituting this relation into Eq.(25), one finds the HDE equation of state parameter 



w A (t) 



-2s d 



kK, 2 6 n 2 3 



1 + 2\f Vt A cosy 



3b 2 (l + Q k ) 



(2? 



According to relation y = y/kRh/a, cosy = 1 when k = 0, so in this case il k — and 
therefore, in flat universe, the HDE equation of state is given by 



W A (t) 



L + 2v^I+^ 



(29) 



From Eqs.(6, 27, 28), we have the effective equation of state as 



w 



eff 



2 V) 



A r 2 «*-^«V 3 

or expressed in terms of the redshift z as 

w A JJ (z) = 



- (l + 2\/tt A ~cosyJ 



8^(e^Vo 2 ) 



- + 2 cosy) . 



e -24^ln(l+z) _ ^ K 20 Q 2 3 

On the other hand, the effective equation of state for CDM is given differently by 

eff , , 6 2 fi + a.) 



(30) 



(31) 



(32) 



Now we are in a position to derive two coupled equations whose solutions determine the 
effective equations of state (as in Ref.[24]). Eq.(5) leads to one differential equation for 

-(i + z )^ = -3fi A (i - n A + n k ) (uf - < ff ) + n k n A (i + 3^ 



(33) 



The remaining differential equation for f2 k comes from the derivative of r k in Eq.(15) using 
Eq.(14) as 



-(1 + z 





dz 



eff , ,eff 



fl k (l + fik)(l + 3o;f). (34) 



In order to obtain a solution, we have to solve the above coupled equations numerically 
by considering the initial condition at the present time: ^r\ z=0 > 0, f2 A = 0.73, Q^ =1 = 
0.01 and fi° =0 = 0.0. 
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2.2 Phantom field 



In this subsection we consider a phantom field with a non-minimal coupling, that is, a 
field with an opposite sign in the kinetic term in the Lagrangian [4]. Such models are 
widely used in order to have w\ less than —1. The action of the universe in this case is 



S 



d 4 Xy/—g 



2^ 



l -^a 2 R+ l -g^d,ad v a + X Cm 



(35) 



In this action we have added a phantom field a, which in non-minimally coupled to the 
curvature with coupling parameter £ CT . Lastly, the term Cm represents the matter content 
of the universe and the term \ multiplying it accounts for the interaction. The presence 
of the non-minimal coupling leads to the effective Newton's constant: 



8nG 



eff 



2_2\-! 



(i-e^v) 



(36) 



We shall follow a procedure similar to the one used in the previous subsection in order 
to obtain the equation of state. 

The cosmological equations and the evolution equation for the interacting phantom 
field are given by 



K 2 (p ri 



PA 



1 & 2 + 6^Haa) 



3(1 



a 



a + 3Ha-6Ca[H + 2H z )a = 

Pm + Pa + 3H (p m + pA+p m + Pa) 
and the non-flat HDE can be expressed as 



Q 
-- o 



(37) 

(38) 
(39) 



p A = -(l-£^V)L 



(40) 



where we have used the effective nature of the Newton's constant (36). 

We examine power-law solutions of equations (37)-(39), in the case of a dark-energy 
dominated universe (p m ,p m 1). Thus, we impose: 



a(t) = a t r 
a(t) = a t s 



Using (40) we have that 



PA (t) = - 2 {1-^k 2 a 2 t 2s °)L- 2 



(41) 



(42) 



and as in the previous subsection L is defined as in Eq.(17). By taking the definition of 
Eq.(42) for the HDE density p\(t), and making use of Eqs.(22, 23) one can obtain 



Pa 



-2s a t Va 



£ a K 2 a z t 



2+2s a 



- 2H (l - a/T^cos?/) Pa- 



(43) 



1 - £ a K 2 a H 2s - 

Substitution of the above relation into Eq.(3) and use of the definition Q = TpA, yields 



w A (t) 



t 



-2s„ 



i 2v^a" r 

— I cos y H 

3 3 y 3H 



(44) 



Inserting Eq.(27) into Eq.(44) allows us to obtain the equation of state parameter 



w A (t) 



t~ 2Sa — £o-K 2 cro 2 



1 + 2y fl\ cosy 



3b 2 (l + Q k 



(45) 



Considering the relation y = \fkRh/a, cosy = 1 when k = 0, i.e. ilk = and therefore, 
in flat universe, the HDE equation of state is expressed as 



w A {t) 



t 2s ° - £ ct k 2 <V 3 V ^A, 



(46) 



From Eqs.(6, 27, 45), we arrive at the effective equation of state 

eff _ 8U^W) 1 



W 



-2s a 



1 + 2y f^A cosy 



(47) 



Similarly to what was done in the previous subsection, we can express (47) in terms of 
the redshift z obtaining 



w A JJ {z) 



e 2^ln{l+z) _ £ aK 2 aQ 2 



- (l + 2^/Q^cosy^j . 



(it 



Finally, we insert Eq.(48) into Eqs.(33) and (34) and solve them numerically by considering 
the initial condition at the present time: ^j-| 2=0 > 0, £l A = 0.73,f2° =1 = 0.01 and 
^=o = 0-0. 



2.3 Quintom model 

In this subsection we consider the quintom cosmological scenario [5], where we consider 
simultaneously a canonical and a phantom field, both with non-minimal coupling. As we 
have stated in the introduction, this combined cosmological paradigm has been shown to 
be capable to describe the crossing of the phantom divide w\ = — 1. The action of this 
model is given by 



S 



d 4 Xy/—g 



±R - l -^R - \^R- 



M 



(49) 



2 J T T 

and the presence of the non-minimal coupling leads to the effective Newton's constant 

87rG e// = k 2 [1- K 2 (^(j) 2 + i a a 2 )] _1 . (50) 



The cosmological equations and the evolution equation for the canonical and phantom 
fields in the interacting case are the following 



H 2 



K 2 [pr 



Pa 



¥ 2 



1 ~2 



3[1 



+ 6^ 2 )] 



3H(f) + 6^[H + 2H 



-Q 



(51) 
(52) 
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a 



a + 3H& - 6£ CT (H + 2H 2 ) a 



-Q 



(53) 



Pm + Pa + 3# (p m + p A + p m + p A ) = 
as usual the non-flat HDE is given by 



(54) 



(55) 



after making use of the effective nature of the Newton's constant (50). We examine power- 
law solutions of equations (51)-(54), in the case of a dark-energy dominated universe 
(Pm,Pm *C 1). Thus, we impose: 



a{t) = a t r 

m = u s " 

a(t) = a t s °. 



Substituting into (55) yields 

PA(t) = ^ [1 



(56) 



(57) 



and as in the previous subsections L is defined as in Eq.(17). 

Using the definition of Eq.(57) for the HDE density p A (t), and considering Eqs.(22, 23) 
gives us 



Pa 



-2 ( S(P + s a ) r 1 [k 2 {^<Plt 2s * + ^cr 2 t 2sa )} -2H(l- y^cosy^j p A . (58 



The substitution of the above relation into Eq.(3) and the use of the definition Q = Tp A , 
yields 



w A (t) = 16 



* 2 (ffiSt 2 '*-&# 2 ") 



1 2y^U 

- H cos y , 

3 3 y 3F 



T777 • (59) 



Inserting Eq.(27) into Eq.(59) gives the equation of state parameter 



w A (t) = 16 



* 2 (ggg^ - exgj 

^ 2 (^ 2 ^+^ 2Sct )-1 



1 + 2y^A cosy + 



3& 2 (i + o fc ; 



(60) 



Given the relation y = \^kRh/a, cosy = 1 when k — 0, i.e. £7^ = and therefore, in flat 
universe, the HDE equation of state is given by 



w A (t) = 16 



« 2 

ft 28 * + ^V 2 t 2s °) - 



(61) 



From Eqs.(6, 27, 60), we have the effective equation of state as 



w 



eff 



16 



* 2 (m^ - ey t 2s °) 
K 2 - 



- (l + 2^/tt A ~cosy^ . 



(62) 
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that can be expressed in terms of the redshift z as 



w 



eff. 



z) = 16 



« 2 {e^le 2 ^ ln{1+z) - ga 2 e- 24 ^+z)) 
k 2 (^0oe 24 ^ /n(1+2) + ^a 2 e- 24 «-^( 1 + z )) - 



- (1 + 2^ cosy). (63) 



Finally, as in previous subsections, we insert Eq.(63) into Eqs.(33) and (34) and solve 
them numerically by considering the initial condition at the present time: 
0, tt° A = 0.73, fig =1 = 0.01 and fi£ =0 = 0.0. 



(j^A I 
dz 



=o > 



3 Cosmological implications 

In the previous subsections we have obtained the equation of state parameter of dark 
energy, w e ^ (z), in terms of the coupling parameters £ CT and the amplitudes 4>q, gq. In 
the present section we investigate the cosmological implications for each case. 

3.1 Canonical field 

In the case of an interacting canonical field, non- minimally coupled to gravity, iw^ (z) 
is given by the relation (31). In Fig.l we depict w e /^ (z) for two different values of the 
coupling and for three different values of the combination k 2 0q- Note that the physical 
requirement of an expanding universe results in an upper limit for £0, namely ^ < 1/6 
(see [26]). 

As we observe, the value of w A (z) at z — 0, that is, its current value w^f, decreases as 
£0 increases, while its dependence on ft 2 ^ 2 , is non-monotonic. However, in this interacting 
canonical field case 10^0 * s always greater than —1, independently of the values of £^ and 
k 2 0q. This was expected since this case is known to be insufficient to describe the crossing 
of the phantom divide ?% = — ^ fr° m above [28]. 

Secondly, we can see that for k 2 0q °f the order of 1 or below, we obtain a divergence of 
w e j[ . This behavior is a clear prediction of relation (31), since it possesses a singularity 
at 

^ = -l + (^ 2 2 )"^. (64) 

Therefore, the combinations of £^ and k 2 0q that satisfy the equation giving a positive z s 
must be excluded. 

Finally, we mention that the effect of non-flat universe is negligible as the curves for 
k = +1,0 appear superimposed. 

3.2 Phantom field 

In the case of an interacting phantom field, non-minimally coupled to gravity, w e /^ \z) is 
given by relation (48). In Fig. 2 we depict w e /^(z) for two different values of the coupling 
£ CT and for three different values of the combination k 2 <t 2 . Note that in this case the 
physical requirement of an expanding universe, results in an upper limit for namely 
-1/6 < & (see [26]). 

As we can see, the value of w e / Q (z) is now a non-monotonic function of £ CT and k 2 g^. 
Furthermore, we observe that for some particular combinations of £ CT and K (Jq j clS £L 
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consequence of the singularity of (48), there is a divergence of (z) at 

z s = -l+(^K 2 a 2 y^ . (65) 

Thus, the combinations of £ CT and k 2 ct 2 that satisfy this transcendental equation giving a 
positive z s must be excluded. 

In the case at hand we can see that w e /J(z) is always greater than —1, independently 
of the values of £ CT and n 2 a 2 . This is expected as we cannot have w e /J(z) < —1 for a 
phantom field in interacting HDE (see [32]). Furthermore, the effect of non-flat universe 
is negligible as the curves for k = +1,0 appear superimposed. Therefore, the non-flat 
universe cannot induce the phantom phase even if one includes a non-minimal coupling 
in the interacting HDE framework. 

3.3 Quintom model 

In the case of the combined quintom model, that is, when both the canonical and phantom 
fields are considered to be non-minimally coupled to gravity simultaneously, w e /^(z) is 
given by relation (63). In Fig.3 we depict w e /^ (z) for two different values of the coupling 

and for three different combinations of k 2 0q an d k 2 (t 2 - Note that in this case the 
physical requirement of an expanding universe, results in an upper limit for namely 

< 1/6 (see [26]). The value of w^f(z) is a monotonic function of As in the previous 
cases, for some particular combinations of n 2 (pl and n 2 a 2 , as a consequence of (63), 
there is a singularity of w e /^(z) at a specific z s . The form of the denominator of (63) 
does not allow for an explicit expression of z s , but a numerical investigation provides the 
specific excluded parameter values. 

As we can observe in Fig. 3, w e / ^(z) is greater than —1, independently of the values of 
k 2 0o an d k 2 <j 2 . Once again, the effect of non-flat universe is negligible as the curves 
for k = +1,0 appear superimposed. It turns out that not even the quintom scenario with 
non-minimal interacting HDE can describe the phantom regime. 

4 Conclusions 

Currently scalar fields play crucial roles in modern cosmology. In the inflationary sce- 
nario they generate an exponential rate of evolution of the universe as well as density 
fluctuations due to the vacuum energy. It seems that the presence of a non-minimal 
coupling (NMC) between scalar field and gravity is also necessary. There are many the- 
oretical evidences that suggest the incorporation of an explicit NMC between the scalar 
field and gravity in the action [33]. The NMC arises from quantum corrections and it is 
required also by the renormalization of the corresponding field theory. Amazingly, it has 
been proven that the phantom divide line crossing of dark energy described by a single 
minimally coupled scalar field with general Lagrangian is even unstable with respect to 
the cosmological perturbations realized on the trajectories of the zero measure [34]. This 
fact has motivated a lot of attempts to realize the crossing of the phantom divide line by 
a equation of state parameter of the scalar field used as dark energy candidate in more 
complicated frameworks. 

Studying the interaction between dark energy and ordinary matter will open the pos- 
sibility of detecting dark energy. It should be pointed out that evidence was recently 
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provided by the Abell Cluster A586 in support of the interaction between dark energy 
and dark matter [35]. However, despite the fact that numerous works have been carried 
out, at present there are no stringent observational bounds on the strength of this interac- 
tion [36]. This weakness to set stringent (observational or theoretical) constraints on the 
strength of the coupling between dark energy and dark matter stems from our unaware- 
ness of the nature and origin of the dark components of the Universe. It is therefore more 
than obvious that further work is needed in this direction. Due to this, we have extended 
our work in [26] to the interacting case in this paper. As a result, in the present paper we 
have investigated canonical, phantom and quintom models, with the various fields being 
non-minimally coupled to gravity, in the framework of the interacting HDE in non-flat 
universe. For this case, the characteristic length is no more the radius of the event hori- 
zon (Re) but the event horizon radius as measured from the sphere of the horizon (L). 
In each case we have extracted , that is, the dark energy effective equation of state 
parameter, as a function of the redshift and used as parameters the couplings and the 
amplitudes of the fields. Finally, we have analyzed it in order to obtain its cosmological 
implications. 
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Figure 1: vs z in the interacting canonical field case, for £^ = 1/7, ^ = 1/9, 

b 2 = 0.01 and k = +1,0, where in each case the combination k 2 0o is taken equal to 10, 1, 
0.1 respectively. The curves for k — +1,0 appear superimposed showing that the effect of 
the non-flat universe is negligible. The divergence of w e / is a direct consequence of the 
singularity of (31 ), and thus the corresponding combinations of and k 2 4>1 must be excluded. 
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Figure 2: itf^ z in the interacting 'phantom field case, for £ CT = 1/4, £ CT = 1/7, 
6 2 = 0.01 and k = +1,0, where in each case the combination k 2 o~q is taken equal to 3, 1, 
0.1 respectively. The curves for k — +1,0 appear superimposed showing that the effect of 
the non-flat universe is negligible. The divergence of w e / is a direct consequence of the 
singularity of (48), and thus the corresponding combinations of and k 2 o~q must be excluded. 
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Figure 3: h>a vs z in the combined interacting quintom scenario, for^ = l/6 + , £</, = 1/8, 
6 2 = 0.01 and k = +1,0, where in each case the combinations k 2 0q an d k2 °~o are shown 
in the insets. The curves for k = +1, appear superimposed showing that the effect of 
the non-flat universe is negligible. The divergence of ^ s a direct consequence of the 
singularity of (63), and thus the corresponding combinations 0/ £</,, k 2 0q and k 2 cj 2 must 
be excluded. 
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